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FOREWORD 


This report was prepared by the Aerospace Sciences Laboratory of 
the Lockheed Missiles and Space Company, Sunnyvale, California. It pre- 
sents the documentation for a digital computer program known as PRESTO II. 
This program was developed by Lockheed for the Langley Research Center 
under NASA Contract NAS 1-5255* A FORTRAN source deck listing included 
with the master copy of this report completes the program documentation. 
PRESTO TT makes use of computer subroutines that were developed for 
PRESTO (Program for Rapid Earth- to- Space Trajectory Optimization) under 
NASA Contract NAS 1-2678. r rhis report does not include any information 
on those subroutines that were described in the PRESTO final report. 

R. C. Rosenbaum was responsible for the work reported here. Computer pro- 
gramming was performed by Mss Zoe Taulbee. 


Results of the program carried out under Contract No. NAS I -2678 
have been published as NASA CR-I 58 entitled "Program for Rapid Earth-to- 
Space Trajectory Optimization" by Robert E. Willwerth, Jr., Richard C. 
Rosenbaum, and Wong Chuck. 
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ABSTRACT 


A high-speed computer program for trajectory optimization has been 
developed. Program speed is obtained through the use of an approximate 
analytic solution to the three-dimensional equations of motion on a 
"cylindrical" Sarth. The thrust attitude control history is assumed to 
follow the "linear- tangent" law. The method of steepest descent is used 
to determine the control parameters that maximize payload while simul- 
taneously satisfying terminal constraints. 
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1. INTRODUCTION 


♦ 


The development of the high-speed digital computer has made it possible 
to apply sophisticated numerical techniques to the problem of trajectory 
optimization. The steepest-descent method, in particular, has been used 
quite successfully to evaluate the performance of rocket boosters. For 
preliminary design studies, however, where a great many configurations 
must be evaluated, conventional numerical methods may require an excessive 
amount of computer time. Most of this time is taken up by repeated inte- 
gration of numerous differential equations. 

In PRESTO II, most of the numerical integration is eliminated through 
the use of an approximate analytical solution to the equations of motion. 

The pitch and yaw components of the thrust attitude history for the upper 
stages are made up of one or two segments which follow the linear-tangent 
rule, i.e., the tangent of the attitude angle is equal to a + bt. The first 
stage is assumed to fly at zero angle of attack. An iterative procedure, 
employing the steepest-descent method, is used to find the values of the 
controls which meet terminal constraints and maximize payload. In addition 
to the thrust attitude history, the controls that can be optimized are the 
initial flight path angle and azimuth, the launch time, the length of a 
coast between two powered stages, and the cutoff and ignition times for a 
restartable final stage. 

It is known that the linear- tangent law provides an optimal solution to 
the ascent problem if one makes a "flat- Barth” assumption. Furthermore, a 
closed-form solution to the equations of motion is available for the M flat- 
Earth" case. In PRESTO II, the Earth is represented by a cylinder. The 
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equations of motion in the pitch plane include the centrifugal and Coriolis 
accelerations. Lateral motion is assumed to be uncoupled from the pitch- 
plane motion. The linear- tangent law is applied to a thrust attitude angle 
that is measured with respect to the local horizontal rather than a fixed 
direction. 

The trajectory of the first stage of the booster vehicle is integrated 
numerically. The first stage trajectory is a function only of the initial 
flight path angle and launch weight because of the zero angle of attack 
assumption. A table of trajectory variables at the end of the first stage 
as a function of flight path angle is formed. This table is then used to 
represent the first stage during the optimization process. 

After an optimum trajectory has been obtained on PRESTO II, the optimum 
controls are used to run a numerically integrated trajectory on PRESTO. 

Terminal conditions will, in general, not be satisfied because of the approxi- 
mations in the closed- form solution. A backward guidance run and a final 
guided run are then made on PRESTO to produce a numerically integrated tra- 
jectory which meets terminal conditions and provides a near-optimum payload. 

A description of the program, together with flow diagrams for the larger 
subroutines, is presented in Section 2. The nomenclature introduced in the 
PRESTO II subroutines is defined in Section 3» Section 4 serves as a manual 
for operating the program. The optimization equations are derived in Section 5 
and Section 6 provides a derivation of the approximate solution to the equa- 
tions of motion. 

PRESTO II makes use of the nomenclature and logical arrangement of 
PRESTO to as great a degree as possible. It is therefore suggested that the 
reader have the PRESTO final report available for reference. 
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2. PROGRAM DESCRIPTION 


PRESTO II is written in the Fortran IV language. The program is composed 
of two distinct parts. One part is the DELTA version of the original PRESTO 
program with a few subroutines altered. Ihe other part is composed of the 
new subroutines associated with PRESTO II. In this section, the two parts 
of the program will be referred to as PRESTO and PRESTO II. 

PRESTO II is comprised of the following subroutines: 


MAIN 
STAGE 1 
PERTRB 
TRAJ 2 
NSTOP 2 
COAST 2 
CONC 
MEQ 2 


CONVRT 
VALUAT 
LOOK 2 
OPCNL 2 
RUNGEK 
NTSGRT 
LIFT 
EQUATN 


The following subroutines are taken from the PRESTO program: 


PRESTO (MAIN of PRESTO) MISCON * 

TRAJ SYMVRT * 

PCAL LOOK * 

OPCNL INER * 

MEQ ATMOS * 

INSTOP UINEPH * 

ICS PLANEP * 

DEQ REIN * 

COAST OPTION * 

RKAD WIC * 


The function of each of the PRESTO II subroutines is discussed here. 
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Subroutine Functions 


MAIN 

The operations performed in the PRESTO II MAIN program are similar to 
those performed in the PRESTO MAIN program. Additional calculations are 
made to set up the initial values of the control and vehicle parameters. The 
closed-form lift-off calculation has been removed to a separate subroutine. 

At the end of the PRESTO II iterations, subroutine PRESTO is called to begin 
the numerical integration of the trajectory. 

STAGE 1 

The trajectory for the first stage of a vehicle launched from the I3arth's 
surface is computed here. A table of trajectory variables at the end of the 
first stage, as a function of the initial flight path angle, is stored in 
the two dimensional array, STABLE. At the beginning of each case, two entries 
are made in GTABLE. The first is for the initial flight path angle, Y^ , and 
the second is for Yj plus DELGAM. The slope of the trajectory variables is 
then computed and stored in QGTABL. Thereafter, new entries in GTABLE and 
QGTABL are made only when the initial flight path angle falls outside of 
the region already included in the table. There is room for twenty-five 
entries on each side of the initial table entry. 

If the fixed-final stage option is being used, so that the launch weight 
is adjustable, the influence of launch weight changes must also be determined. 
Initially, this is done by perturbing the launch mass by .001 times the 
initial mass. The change in the first stage burnout conditions due to this 
mass change is stored in the OKMS array. There is an entry in OKMS for every 
entry in GTABLE. Thus, every initial flight path angle has a mass slope 
associated with it. Furthermore, for every table entry, the launch mass is 
stored in the array GMASS. When a table lookup is made, the current launch 
mass is compared with the GMASS associated with the current flight path 
angle. This mass difference is multiplied by the current value of the mass 
slope, OKM 1, to obtain the correction in first stage burnout conditions due 
to the change in launch weight. A new mass slope is computed after NEWKM 1 
entries have been made in GTABLE. NEWKM 1 is input by the user. 
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The same sort of correction is made if the launch azimuth is adjustable. 
The current value of the azimuth slope is stored in the second half of OKM 1. 

A new slope is computed whenever the initial azimuth changes by more than 
DELPSI. DELPSI is also input by the user. 

The trajectory variables at the end of the first stage are placed in FI 
and are returned to the calling program by means of an argument list. 

A flow chart for STAGE 1 is at the end of this section. 

PERTRB 

The sensitivities of the terminal constraints and stopping parameters to 
changes in the control variables are computed here. PERTRB replaces the back- 
ward integration performed in PRESTO. Each of the controls is perturbed by 
an amount, DKLCON, set by the user. A trajectory is then computed. The 
perturbed terminal conditions are compared with the nominal terminal conditions 
to determine the influence of the control change. 

A flow chart for PERTRB is at the end of this section. 

TRAJ 2 

TRAJ 2 performs the same operations as does TRAJ in PRESTO. TRAJ 2 is 
considerably smaller than TRAJ because of the absence of a backward trajectory. 
The subroutine has also been shortened through the use of a common area for 
determining the trajectory during each stage. 

A flow chart for TRAJ 2 is at the end of this section. 

N3T0P 2 

This subroutine interpolates to the correct stopping parameter. 

COAST 2 

Motion across a coast stage is determined here. The equations used are 
from the COAST routine in PRESTO. 

OONC 

The control variables to be used on the next trajectory are computed. 
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MEQ 2 

The matrices required to compute the control changes in CONC are com- 
puted. See Section 5 for a derivation of the equations used. 

CONVRT 

The trajectory variables associated with the cylindrical coordinate 
system are converted to the geographic coordinates used in PRSSTO. The equa- 
tions used are given in Section 6. 

VALUAT 

The derivatives used in the numerical integra+.i on of the first stage 
are computed. 

LOOK 2 

This routine is similar to LOOK. It provides the linear interpolation 
used in STAC? 1. 

OPCNL 2 

The trajectory variables at the beginning and end of each stage are 
output. 

RUNOI ffC 

This subroutine provides the numerical integration formulas used for the 
first stage. It is similar to RKAD. 

NTKGRT 

The integration of the first stage is controlled. 

LIFT 

The closed-form lift-off calculation is performed. 

EQUATN 

The approximate analytic solution to the equations of motion is computed. 
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The results are stored in the array, UVH, and are returned to the calling 
program through an argument list. 

The following PRESTO subroutines have been changed. 

PRESTO 

PRESTO is the subroutine name for the MAIN routine of the PRESTO program. 
The initial setup operations now appear in the PRESTO II MAIN and have been 
removed from subroutine PRESTO. When the program is under the control of 
PRESTO, it runs a nominal trajectory, a backward guidance trajectory, and a 
final guided run. Control is then returned to the MAIN program. 

TRAJ 

A change has been made so that an intermediate constraint (circular park 
orbit or radius of perigee) will be applied on the final guided run if it was 
used in PRESTO II. 

PCAL 

The thrust attitude angles are computed from the PRESTO II angles on 
the nominal trajectory. This computation replaces the table look-up in PRESTO. 
The equations used here are derived in Section 6. 

REIN and WIC 

Changes associated with the PRESTO II data requirements have been made. 

Use of PRESTO II Without PRESTO 

If a numerically integrated trajectory is not required, the size of the 
program can be greatly reduced by eliminating many of the PRESTO subroutines. 
PRESTO II can run by itself using all of the PRESTO II subroutines plus those 
PRESTO subroutines identified by an asterisk on page 2-1. A dummy subroutine 
PRESTO is also required. 
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Use of COMMON Variables 


'Jbe variables stored in COMMON have been separated into three groups. 

One group, known as the common COMMON, is used in all subroutines that employ 
COMMON. A second group of variables which are used only in PRESTO II sub- 
routines appears at the end of common COMMON in the PRESTO II subroutines. 

A third group of variables which are used only in the PRESTO subroutines 
appears at the end of common COMMON in the PRESTO subroutines. 
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3- DEFINITIONS OF PRESTO II VARIABLES 


The variables introduced in the PRESTO II subroutines are defined here. 


CA 

CB 

CE 

CF 

CON 


CNTRL 1 
CNTRL 2 


Pitch plane steady- state control angle for current stage 
Pitch plane slope for current stage 

yaw plane steady-state control angle for current stage 

Yaw plane slope for current stage 

Vector of current values of adjustable parameters 

3* Length of coast after first burn in stage 4 

4. Length of first burn in stage 4 

5. Length of coast after stage 3 

6. Length of coast after stage 2 

7* Length of coast after stage 1 

8. Launch time of day 

9. Initial azimuth 

10. Initial flight path angle 

11. Initial value of pitch plane control variable for segment one (a^) 

12. Slope of pitch plane control variable for segment one (b^) 

13. Initial value of yaw plane control variable for segment one (e^) 

14. Slope of yaw plane control variable for segment one (f^) 

15. Initial value of pitch plane control variable for segment two (a 2 ) 

16. Slope of pitch plane control variable for segment two (b^) 

17* Initial value of yaw plane control variable for segment two (e^) 

18. slope of yaw plane control variable for segment two (fg) 

Data block 40 

Data block 4l 


DSLCON Vector of magnitudes of perturbations of adjustable parameters 
DELGAM Bata block 43 
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DKIMAS 

DKLPSI 

DELTA 

DTR 

EQC 

EQT 

EXTRA 


PI 

GAMONE 

GMASS 

GPS I 

GTABLE 

HAV 
HAVI 
IS EG 
J6 


Magnitude of launch mass perturbation 
Data block 43 
Data block 42 

Conversion factor for converting degrees to radians 

Vector of powered stages exhaust velocities 

Vector of thrusts for powered stages 

Vector of miscellaneous variables 

1 Stored value of ST3(l) preceding circular park orbit 
2-7 Stored values of unperturbed stage 1 burnout conditions 

used to determine sensitivity of stage 1 burnout conditions 
to initial azimuth 

8 Stored value of GAMI at start of the circular park orbit 

9 Stored value of X(3) at start of the circular park orbit 

10 Stored value of X(4) at the end of the trajectory 

Vector established in STAGE 1 routine with stage 1 burnout values 
of u, v, h, T, and X, and returned to calling routine 

Initial flight path single 

Launch mass vector associated with GTABLE row entries (for initial 
flight path angle y ) of integrated first stage final trajectory 
variables u, v, h, t, and X (in STAGE 1 routine) 

Initial azimuth vector associated with GTABLE row entries (for 
initial flight path angle Yj) of integrated first stage final 
trajectory variables u, v, h, T, and X (in STAGE 1 routine) 

Pinal trajectory variables u, v, h, t, ^ , and X for values of 
initial flight path angle Vj* A table matrix established and used 
by STAGE 1 routine 

Vector of average altitudes of powered stages 
Data block 45 
Data block 44 

Used in MBQ2 to indicate whether optimization indicator shall be computed 
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JSTABL 


Code of adjustable parameter being processed in PERTRB routine 


KASHI Used in STACK 1 routine to indicate whether initialization of GTABLE 

is needed 

KNSTRN Code of terminal constraint being processed in PERTRB routine 

NADJST Sequence number of adjustable parameter being processed in PERTRB 
routine 

NCNSTR Sequence number of terminal constraint being processed in PERTRB routine 
NEWKM1 Data block 43 

HUMS EG A vector indicating for each stage its relation to control segments: 

1 for being included in control segment one; 2 for being in control 
segment two; and 0 for not being included in a control segment 

IiLl First powered stage in segment one (KK4) 

LI2 Last powered stage in segment one (KK4) 

LL3 First powered stage in segment two (KK4) 

LlA Current segment number 

LL5 Number of segments 

Ll£ Last stage used in a forward trajectory (KKl) 

LL7 First powered stage in segment one (KKl) 

LL8 First powered stage in segment two (KKl) 

LL9 Adjustable coast stage (KKl) 

LL10 First burn in fourth stage (KKl) 

LT.11 Number of angle controls plus other adjustable parameters 

LLI2 MISCON calling key 0 - do not call MISCON 

1 - call MISCON 

LL13 Adjustable coast stage (KK2) 

T.L14 Last powered stage in trajectory (KK4) 

I.L15 Segment indicator of adjustable controls being perturbed in PERTRB 

T.L16 Indicates presence of stage 8 0 - stage 8 not used 

1 - stage 8 used 
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LL17 

i.li8 


LL 19 

LT20 

U2.L 

LTJ22 

1.12 3 

LI24 

LI25 

LT26 

LT27 

U28 

LL29 

LL30 

0KM1 

OKMS 

OKS 

OKZ 


Available 

Trajectory parameter printout key 1 - initial stage 1 values 

2 - final stage 1 values 

3 - initial powered stages 

4 - final powered stages 

5 - initial coast stages 

6 - final coast stages 

Weight (OMI i ) subscript for printout (0PCNI2) 

Nominal or perturbed trajectory 0 - nominal trajectory 

1 - perturbed trajectory 

Current value for the powered stage immediately following the last 
coast stage. If there is no coast stage, it is equal to LL7- 

Assigned go to integer variable for LI21 

Available 

Available 

Available 

Available 

Available 

Available 

Available 

Available 

First six elements are sensitivities of stage 1 burnout conditions 
to launch mass. Other six elements are sensitivities of stage 1 burn- 
out conditions to initial azimuth. 

Stored values of first six elements of 0KM1 

Sensitivity matrix of terminal constraints to adjustable parameters 
for a fixed value of the stopping parameter 

Sensitivity matrix of terminal constraints to adjustable parameters 
for a fixed value of the final stage mass ratio 


OLS 

OUZ, 

OMI 

OMZ 

ONZ 

PCON 

PGAMA1 

POMI 

PPAY 

PSA 

PSB 

PSE 

PSF 

PS IONS 
PSIPR 

PSTGTM 


Sensitivity of terminal constraints to launch mass for a fixed value 
of the stopping parameter 

Sensitivity of terminal constraints to launch mass for a fixed value 
of the final stage mass ratio 

Vector of each powered stage's initial and final masses 

Sensitivity of stopping parameter to adjustable parameters for a 
fixed value of the final stage mass ratio 

Sensitivity of stopping parameter to launch mass for a fixed value 
of the final stage mass ratio 

Vector of values of achieved adjustable parameters which are perturbed 
by the vector DELCON in the PERTRB routine 

Value of initial flight path angle used by PERTRB routine for compu- 
tations of perturbed trajectories 

Values of initial and final masses of powered stages used by PERTRB 
routine for computations of perturbed trajectories 

Sensitivities relating final weight to errors in terminal constraints 

Vector of pitch plane steady- state control angles for control segments 
one and two used by PERTRB routine for computations of perturbed 
trajectories 

Vector of pitch plane slopes for control segments one and two used 
by PERTRB routine for computations of perturbed trajectories 

Vector of yaw plane steady- state control singles for control segments 
one and two used by PERTRB routine for computations of perturbed 
trajectories 

Vector of yaw plane slopes for control segments one and two used by 
PERTRB routine for computations of perturbed trajectories 

Initial azimuth 

Vector of derivatives of terminal constraints with respect to last 
stage mass ratio 

Vector of stage time durations for all stages used by PERTRB routine 
for computations of perturbed trajectories 
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PUVH 

PZETA 

QGTABL 

RHAV 

SA 

SB 

SE 

SEGTIM 

SF 

SF1 

SHAV 

SOMI 

SPR 

SST3 

STGTIM 

SWFUEL 

TABLET, 

TCOAST 

TCTIM 

TIM 

UVH 


Matrix of trajectory variables u, v, h, t, ^ , X, P, v, and A for 
initial and final values of each stage used by PERTRB routine for 
computation of perturbed trajectories 

Vector of powered stages mass ratios used by PERTRB routine for com- 
putation of perturbed trajectories 

Matrix of slopes associated with GTABLE matrix entries in STAGE 1 
Vector of average radii of powered stages 

Vector of pitch plane steady-state control angles for control segments 
one and two 

Vector of pitch plane sbpes for control segments one and two 

Vector of yaw plane steady- state control angles for control segments 
one and two 

Vector with values of control segment time duration up to that stage 
for each stage included in a control segment, with a zero stored in 
stages not included in a control segment 

Vector of yaw plane slopes for control segments one and two 

Vector returned by STAGE 1 routine with stage one burnout values of 
u, v, h, t, ijr and X 

Vector of saved values of HAV vector 

Vector of saved values of OMI vector 

Derivative of stopping parameter with respect to mass ratio 

Vector of saved values of ST3 vector 

Vector of stage time duration for all stages 

Correct fuel weight in the final stage 

Data block 43 

Time duration of a coast stage 

Total time in each control segment initially 

Current value of control segment time 

Matrix of trajectory variables u, v, h, t, Jr , p, w and A for initial 
and final values of each stage 
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ZETA 

ZINV 

ZKA 

ZKAL 


Vector of mass ratios for each powered stage 
Inverse of weighting function matrix 

Matrix relating changes in terminal constraints to changes in 
control parameters 

Vector relating changes in terminal constraints to changes in launch 
mass 
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4. PROGRAM OPERATION 


DATA INPUT FORMAT 


Data Blocks 


Data required for PRESTO execution are grouped and input in data blocks which 
are identified by number. Each block contains a common type of information 
such as the case title or a stage thrust table, and utilizes a specified 
FORTRAN format for the entire data block. The contents of each data block 
are discussed and then summarized on the next several pages. 

Header Cards 


Input data are punched on cards which are placed behind the program binary 
deck. Cards for each data block are preceded and identified by a "header" 
card. The format of the header cards is 413, where the first field is the 
data block number and the second and third fields give the (inclusive) 
locations within the data block that the subsequent data cards are to be 
placed. Thus, if the user wished to change constants 2, 3 and 4 in data 

block 10, the header card would be punched 10 2 4 , and 

the first three fields of the next card would contain the three constants* 
The fourth field of the header card is used to identify the vehicle stage 
number for the data, where necessary. The stage designation should be 
included only for the data blocks as specified on the following pages. 

Card Sequencing 

Cards within each data block must follow sequentially, but data blocks may 
be arranged in any order. A blank card must follow the last data card. 

Successive Cases 


Successive cases can be run on PRESTO with a minimum of additional data 
input. Only the changes In data from the preceding case must be input. A 
blank card must follow the data for each case. 

999 Card 

A card with 999 punched in the first three columns must end the data deck. 
It is to be placed behind the blank card which ends the data for the final 
case. When the READ routine encounters the 999 card, a normal stop is 
indicated to the machine operator, so that the job can be terminated. 
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DISCUSSION OF THE DATA BLOCKS 


The following discussion will emphasize only the changes in the input 

with respect to the PRESTO input. The summary definition of the data blocks 

will be complete. 

Data Block 4 

(4) For planar missions, optimize theta only. For non-planar missions, 
optimize both theta and chi. 

(ll) The q-alpha constraint is no longer available. If the first stage is 
being used, the zero-alpha constraint is automatically applied during 
PRESTO II iterations and must, therefore, be specified in the input so 
that the final PRESTO runs will be correct. Positions 5 and 6 in data 
block 2 2 must correspond to the burn time for the first stage. 

(18) This option indicates whether or not there is a coast stage between 
the two control segments. 

(19) If the accuracy of a numerically integrated trajectory is not required, 
the program will terminate at the end of the last PRESTO II iteration. 

(20) Checkout matrices and the integration of the first stage can be output 
with this option. 

Data Block 5 

Stage 5 is no longer available. Stage 1 is always the numerically integrated 

stage. For orbital launches, stages 2, 3 or 4 may be used as the first stage 

Data Block 6 


■There is only one adjustable burn in stage 4 and only one adjustable coast 
can be specified. Parameter codes 1 and 2 are not available. 

Dat a Block 8 

Allow no more terminal constraints (including mass) than the number of adjust 
able controls. 

Data Bl ock 24 

The nominal value of coast time is specified in addition to the coast angle. 
The program will iterate to the input time on the nominal trajectory. 'ttie 
input angle is an estimate of the coast angle corresponding to the desired 
time. 


Data Block 27 


The recommended weighting functions are the following: 


Name 

Code 

Weighting Function 

Coast 

3,5, 6, 7 


Burn 

4 

u>i 

Launch Time 

8 

10'? 

Azimuth 

9 

10“ 1 

Flight Path Angle 

10 

500. 

Initial theta in first segment 

11 

1 . 

Slope of theta in first segment 

12 

1000. 

Initial chi in first segment 

13 

.1 

Slope of chi in first segment 

14 

1000. 

Initial theta in second segment 

15 

.1 

Slope of theta in second segment 

16 

10 to 1000 

Initial chi in second segment 

17 

.1 

Slope of chi in second segment 

18 

10 to 1000 

Data Block 30 



Thrust must be constant except for the 

first stage. 



Data Blocks 31, 32, 33 


The only aerodynamic coefficient to be used is the drag coefficient for stage 1* 


Data Block 4o 


The first and last values of theta for each control segment are specified. If 
there is no coast between the two control segments, theta must be continuous, 
i.e,, the last value of the first segment must equal the first value of the 
second segment. 

Data Block 4l 


Same as 40 for chi. 
Data Block 4-2 


The recommended value for the perturbation in the average value of a control 
angle is .01 radians. For a slope, the recommended value is .0001 radians/sec. 
If a stage is quite short, a larger perturbation should be used for the slope. 
The following values are suggested for the other perturbations: 

Coast time 5 sec 

Burn time 1 sec 

Azimuth .01 rad 

Flight path angle .001 rad 
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Data Block 43 


DELGAM is the separation between entries in the table for stage 1. Provision 
has been made for 25 entries on each side of the initial value. A value of 
.1 degree is stiff icient for vehicles with an initial flight path angle ranging 
from 87 to 89 degrees. A larger value for DELGAM can be used if the initial 
flight path angle is lower. 

NEWKM1 is the number of gamma table entries in between evaluation of a new 
sensitivity for launch weight perturbations. 3 is a satisfactory value for a 
DELGAM of .1. 

DELPSI is the change in azimuth required before a new sensitivity to initial 
azimuth perturbations is computed. Ten degrees is the suggested value. 

Data Block 44 


This data block determines the number of control segments to be used and the 
stages contained in each. The first number is the KK4 corresponding to the 
first powered stage in segment one. This is followed by the KK4 's of the 
remaining stages in the first segment, including coasts. Six is input at the 
end of the first segment. If there is only one segment, a 7 follows the 6 . 

If there are two segments, the KK4 *s of the stages in the second segment follow 
the 6 and a 7 terminates the data block. 

Example: If the stage sequence in data block 5 is 

16273484 

and two segments are to be used with stage 8 in between 
the two segments, data block 44 should contain 

293^657 

Data Block b5 

An expected average altitude must be input for each stage. 


I 


DATA INPUT 


DATA BLOCK NUMBER 

TITLE 

FORMAT 

COMMENTS 

1 

Main Heading 

12A6 


2 

Case Heading 

12A6 


3 

Date 

2A6 


k 

Options 

2^11 

IP(n) - See list of available options 

5 

Stage Sequence 

1211 

ISTGE(n) 

Powered flight stages 1 through 4 
Closed-form coast stages 6 through 9 
* End sequence with zero 

6 

Adjustable 

Parameters 

1113 

IBl(n) 

IBl(l) = The number of adjustable para- 
meters to be optimized 
IBl(2) * The adjustable parameter codes 
IBl(ll) non “^ ecreasin 8 order 


3 

4 

5 

6 

7 

8 

9 

10 


1 

2 

3 

4 

5 

6 


** Adjustable Parameter Codes ** 


Length of coast after first burn in stage 4 

Length of first burn in stage 4 

Length of coast after stage 3 

length of coast after stage 2 

Length of coast after stage 1 

Launch time of day 

Initial azimuth 

Initial flight path angle 


Terminal Constraints 713 


IB3(n) - Mission-Dependent 
IB3(l) = Stopping parameter code 
IB3(2) = Form or number of terminal 
constraints 

IB3(3) = List of constraint codes for 
l orbit- injection and lunar 

IB3(7) descent missions 


** stopping Parameter Codes ** 

A. Orbit- Injection Mission 


TOOL (Vj 2 - 2 M./R) 

Available 

Available 

Inertial velocity (V^) 
Available 1 

Velocity 
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DATA INPUT 


B. Lunar Transfer Mission 

1 Transfer time to lunar radial distance 

C. Lunar Descent Mission 
Use only Code 6 (Velocity) 

D. Planetary Transfer Mission 
1 Hyperbolic Excess Velocity (VH) 

** Form or Number of Terminal Constraints ** 

A. Orbit-Injection Mission 

and 

C. Lunar Descent Mission 

The number of terminal constraints (l through 6) 

* Include the terminal mass constraint in this count 

B. Lunar Transfer Mission 

1 Use Ephemeris, do not constrain inclination of transfer orbit 

2 Use Ephemeris, do constrain inclination 

D. Planetary Transfer Mission 

1 Do not use Ephemeris (input VH vector) 

2 Mars 

3 Venus 

** Constraint Codes for Orbit- Injection and Lunar Descent Missions ** 

1 Perigee radius 

2 Orbit inclination 

3 Longitude of ascending node 

4 Argument of perigee 

5 

6 Altitude 

7 Inertial flight path angle 

8 Longitude 

9 Inertial azimuth 

10 latitude 

11 Available 

12 Available 

13 Available 

14 Radius 
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DATA INPUT 


DATA BLOCK NUMBER TITLE FORMAT 


COMMENTS 


9 


Output Frequency 2413 


IB4(n) 

Frequency defined as: 

(number of computed points) 
(number of output points) 

IB4(l) - Initial trajectory 

IB4(3) - Final 

IB4(4) - Stage code 

IB4(5) - Initial trajectory 

IB4(7) - Final 

IB4(8) - Stage code 

IB4(9) - Initial trajectory 

IB4(ll)- Final 
IB4(l2)- Stage code 
IB4(13)- Initial trajectory 


IB4(15)~ Final 
IB4(l6)- Stage code 


10 Nominal Constants 11E12.8 CTl(n) 

Nominal values shown 
Input only if changes are desired 
CTl(l) = 7.29211E-5 rad/sec omega 
CT1(2) = 1.407735E16 ft3/sec 2 FMU 
CT1(3) = 20902900. ft RE 
CTl(4) = 1716.4827 gas constant for 
atmosphere subroutine 
CTl(5) = 32.154856 ft/sec 2 measured sea 
level gravity; veight-to-mass 
conversion factor (GO) 

CT1(6) = 2116.2 lb/ft 2 PO 

sea level ambient pres sire 
CT1(7) = 250,000. ft HAERO 

Altitude above which all aero- 
dynamic computations are bypassed 
CTl(8) = 46. Upper limit of lines per 
page 

CT1(9) = (Available) 

CTl(lO)= .02 sec TIMER 

Tlme-epsilon increment to insure 
hitting critical time 
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DATA INPUT 


i 


DATA BLOCK NUMBER TITLE FORMAT 

11 Atmosphere 8E12 . 8 

12 Subroutine 8E12 . 8 

13 Constants 8E12.8 

14 Constants 8E12.8 


15 Nominal Constants 8E12.8 


16 Time Duration 12E12.8 

and Number of 
Integration Steps 
Within Each Stage 


17 Input Data 20E12.8 


19 Terminal 6E12.8 

Constraint 
Magnitudes 


COMMENTS 

CT2(n) Nominal values tabulated 

CT3(n) in report section 

CT4(n) describing atmosphere 

CT5(n) subroutine 

*Input only if changes from nominal 
values are desired 

CT6(n) *Input only for change 
Used in forming partial derivatives 
across closed- form coasts 

CT6(l) = 1.0 Velocity increment 
CT6(2) = .001 Gamma increment 
CT6(3) = 100, Radius increment 
CT6(4) = .001 Tau increment 
CT6(5) = .1 Mass increment 
CT6(6) = .001 Psi increment 
CT6(7) = .001 Lambda increment 
CT6(8) = .001 Coast angle increment 

DAl(n) 

Time in seconds 

DAl(l) Stage 1 time duration 

DAI (2 ) Stage 1 points 

DAI ( 3 ) Stage 2 time 

DAI (4) Stage 2 points 

DAI (5) Stage 3 time 

DAl(6) Stage 3 points 

DAl(7) Stage 4 total time duration 

DAI (8) Stage 4 total points 

DAl(9) Stage 4 first burnout time 

DA2 (n) 

DA2(l) Input should = GO 
DA2 (2) = +1.0 for takeoff 

= - ] . 0 for retro-burns 

DA2(l7) Magnitude for total time 
constraint 

DA4(n) Mission- Dependent 
DA4(l) is the magnitude of the stoppin g 
parameter 

DA4(2) are the desired values of the 
I terminal constraints, listed in 
DA4(6) the same order as in Data Block 
Number 8 

Units are feet, radians, seconds 
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DATA INPUT 


DATA BLOCK NUMBER TITLE FORMAT 


COMMENTS 


** Terminal Constraints Input Guide ** 

A. Orbit-Injection and Lunar Descent Missions 

and 

C. As Described Above 

B. Lunar Transfer Mission 
DA4(l) Transfer time (hours) 

DA4(3) Used internally 

DA4(4) Inclination of transfer orbit (radians) 

D. Planetary Transfer Mission 
DA4(1) Hyperbolic excess velocity 

D a 4 (2 ) Right ascension of hyperbolic asymptote (radians) 
DA4(3) Declination of hyperbolic asymptote (radians) 

DA4(6) Transfer time (days) (DA4(4) and DA4(5) are used 
internally) 

* DA4(l,2,3) are not required input if PLANEP routine is called 


20 Initial Times 2E12.8 DA5(n) 

DA5(l) Time duration for closed- form 
lift-off calculation (sec) 

DA5 (2 ) Launch time in Space Age Date 
(days from January 0.0, i960) 

21 Initial Conditions 6E12.8 DA6(n) 

DA6(1) Altitude (feet) 

DA6(2) Longitude (degrees - from prime 
meridian) 

DA6(3) Latitude (degrees) 

DA6(4) Velocity (feet/second) 

DA6 ( 5 ) Azimuth (degrees) 

DA6(6) Flight path angle (degrees) 


22 

Control Variable 
Constraints 

6E12.8 

DAT(n) 
For a 

= 0 constraint 




DA7(5) 

DA7(6) 

Time to start testing 
Time to release constraint 

23 

Constraints on 
State Variables 
During Boost 

4e12.8 

DAu(n) 

da8(i) 

da8(2) 

da8(3) 

DAo(4) 

Available 

Available 

Minimum altitude for coast 
stages 8 and 9 (feet) 

Aval lab le 
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DATA INPUT 


DATA BLOCK NUMBER TITLE FORMAT 


COMMENTS 


24 


27 


Nominal Range 8E12.8 

Angles and Time for 
Closed-Form Coast 
Stages 


DA9(n) ,n=l,4 Angles in radians of 
central arc 


DA9(I) 

DA9(2) 

DA9(3) 

DA9(4) 


Coast stage 6 
Coast stage 7 
Coast stage 8 
Coast stage 9 


DA9(n)n=5,8 Coast time in seconds 

DA9(5) Coast stage 6 
DA 9(6) Coast stage 7 
DA9(7) Coast stage 8 
DA9(8) Coast stage 9 


Weighting Constants 19E12.8 
for Control Variables 


DA12(n) 

DA12(3) For adjustable parameters 
I identified by code number 

DA12 ( 10 ) as listed for Data Block 6 
DA12(ll) Theta in segment 1 
DA12 ( 12 ) Theta slope in segment 1 
DA12(l3) Chi in segment 1 
DA12(l4) Chi slope in segment 1 
DA12(l5) Theta in segment 2 
DA12(l6) Theta in segment 2 
DA12(17) Chi in segment 2 
DA12(l8) Chi slope in segment 2 


28 


29 


Convergence Data 6E12.8 DA13(n) 

DA13(1) 

DA13(2) 
DA13(3) 
DA13 (k ) 
DA13(5) 

Permitted Values 5E12.8 DAl4(l) 

• 

DAl4(5) 

* Suggest 
one deg 


Initial weight improvement 
to be used if (automatic) inter- 
nal computation fails 
Epsilon- weight for stopping 
attempted mass improvement 
Maximum number of forward tra- 
jectories per case 
Velocity below which Runge 
Kutta integration is always used 
Velocity to start open- loop 
computations 

List in same order as constraints 
are specified and in the same 
units 

10,000 feet for distances and 

ee on angles 


30 


Vacuum Thrust 12E12.8 TT1, TT2, TT3, TT4 

Tables * Stage code required on header cards 

Sequence: Blank, time, thrust, time, 

thrust. !• OfilO 
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DATA INPUT 


DATA BLOCK NUMBER TITLE FORMAT 


COMMENTS 


31 


39 


4o 


4i 


Drag Coefficient 22E12.8 TD1 

Tables (ALPHA = 0) * Stage code required on header cards 

Sequence: Blank, Mach, C D , Mach, C , 

Mach 1.0E10 ° 


Stage Data 


Theta History- 


Chi History 


5E12.8 SGI, SG2, SG3, SG4, each dimensioned 5 
* Stage code required on header cards 

SGx(l) Vacuum I 8 _ (sec) 

SGx(2) Aerodynamic reference area. (ft 2 ) 
SGx(3) Total nozzle exit area (ft'j 
SGx(4) Initial weight (lh) 

SGx(5) Final veigit (lb) 


4E12.8 


* CNTRLl(n) 
CNTRLl(l) 
CNTRL1(2) 
CNTRL1(3) 
CNTRLl(U) 


Units of degrees 
Initial value for first segment 
Final value for first segment 
Initial value for second segment 
Final value for second segment 


4E12.8 


* CNTRI2 (n) 
CNTRI2(1) 
CNTRI2 (2 ) 
CNTRI2(3) 
CNTRI2(4) 


Units of degrees 
Initial value for first segment 
Final value for first segment 
Initial value for second segment 
Final value for second segment 


42 


li.3 

• 


♦When making a change in theta or chi, the 
entire table must be input. 


Perturbations 8E12.8 

Also used in mass 
improvement calcula- 
tion 


DELTA(n) 

DELTA(l) Average value of theta and chi 
in first segment 

DELTA(2) Slope of theta and chi in 
first segment 

DELTA(3) Average value of theta and 
chi in second segment 
DELTA(4) Slope of theta and chi in 
second segment 
DELTA(5) Flight path angle 
DELTA(6) Azimuth 
D£LTA(7) Coast 
DELTA(8) Fourth stage burn 


Stage 1 Controls 


3E12.8 TABLEG(n) 

TABLEC-(l) = DEIGAM in degrees 
TAB LEG (2) = NEWKM1 
TABLBG(3) = DEIPSI in degrees 
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DATA INPUT 


DATA BLOCK NUMBER TITLE FORMAT 


44 Stage Sequence In 1013 

Control Segments 

Stage Number KK4 


2 2 

3 3 

4 - 1st burn 4 

4 - 2nd burn 5 

6 8 

7 9 

8 10 

9 11 

45 Initial Average 5E12.8 


Altitude for Powered 
Stages 


COMMENTS 

ISEG(n) 

KK4 of stages in control segments. 

Place 6 in between segments and 7 at end. 


HAVl(n) 

HAVI(2) Stage 2 
HAVI(3) Stage 3 
HAVI(4) Stage 4 - first burn 
HAVT(5) Stage 4 - second burn 
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DATA INPUT 


** LIST OF AVAILABLE OPTIONS FOR DATA BLOCK NUMBER k ** 


IP(l) Mission 1 

2 

3 

k 

5 

6 


IP(2) 

Axes 

0 



2 

IP(3) 


0 

IP(k) 

Control Variables to 

0 


be optimized 

1 

IP(5) 

Closed- Form Liftoff 

0 


Computation 

1 

IP(6) 

Thrust 

0 

rp(7) 

Mass Computation 

0 

IP(8) 

Aerodynamics 

0 

IP(9) 

Frequency of New Adjoint 
Solution 

1 

IP (10) 

Special Computations and 
Printout 

0 

IP(ll) 

Constraints on Control 

0 


Variables 

2 

IP(12) 

Constraints on State 
Variables 

0 

IP(13) 

Total Time Constraint 

0 



1 

IP(l4) 

Mass -Improvement Procedure 

0 


l 


Orbit- injection 

Inject into lunar transfer from 
Earth liftoff 

Inject into lunar transfer from 
Earth orbit 
Lunar landing 

Inject into planetary transfer from 
Earth liftoff 

Inject into planetary transfer from 
Earth orbit 

Rotating and 3D 
Non-rotating and 3D 


"Diet a only 
Theta and Chi 

Use this computation 
Do not use this computation 

Table look-up - vacuum thrust vs time 

MDOT = V(G 0 * I Bp ) 

Lift = 0. only 

After every successful forward 
trajectory 

None 


None 

ALPHA « 0.0 
None 


Do not constrain total time 
Do constrain total time 

Fixed launch weight 
Fixed final weight 


DATA INPUT 


IP(15) 

IP(i6) 
IP (IT) 

IP (18) 

IP(19) 

IP (20) 


** LIST OP AVAILABLE OPTIONS ** 


Memory Dump 


Output 


0 Do not dump memory 

1 Dump memory at end of job (includes 
floating-point dump of floating-point 
numbers in common) 

1 Non-buffered output 


Output of Input Data 

Coast Between Control 
Segments 

Make Pinal Iterations on 
PRESTO 


0 Do not output data 

1 Output data 

0 No coast between segments 

1 Coast between segments 

0 Do not call PRESTO 

1 Call PRESTO 


Output of Check Matrices 


0 Do not output check matrices 

1 Output check matrices 
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Changes in Program Use as Compared to PRESTO 


There are several changes and limitations in the use of PRESTO II as 
compared to PRESTO. These are listed here. 

1. A cylinder is not a good representation for the Earth if large 
lateral maneuvers are made. The accuracy of PRESTO II will, therefore, fall 
off as the size of the lateral maneuver increases. 

2. The approximation concerning the — term in the equations 

of motion is not good when one stage goes from orbital speed to a high super- 
circular velocity. The program should not be used for interplanetary trans- 
fers requiring a large hyperbolic excess velocity. Program operation is 
satisfactory for lunar transfer missions. 

3. When using the lunar landing mission option, the horizontal 
component of velocity must not be allowed to go negative. This can happen 
if a poor nominal is used or if the stopping parameter is too close to zero. 
Interpolation to the correct stopping condition fails when the velocity goes 
negative. 


4. When the total time constraint option is 

function associated with the final coast stage must be 
_ *8 

e. g. , 10 . 


used, the weighting 
set to a large number, 
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5. DERIVATION AND PROGRAMMING OF OPTIMIZATION EQUATIONS 


The optimization procedure used in PRESTO II is similar to that used in 
PRESTO. The major differences are: 

1. The thrust orientation angles are now discrete rather 
than continuous; therefore no integral terms appear. 

2. The calculation for new controls is made only once at 
the beginning of a trajectory. There is no closed-loop 
calculation as in PRESTO. 

3. The sensitivity of terminal constraints to a change in 
a control is determined by perturbing the control and 
then comparing the perturbed teminal constraint with 
the nominal value. This process replaces the backward 
integration runs in PRESTO. 


The fundamental relation is 

iV - Kdfc + L 


( 1 ) 


where is the terminal constraint vector 

Jfc is the control vector 

d is the perturbation in launch mass 

K is a matrix of sensitivities, or 
influence coefficients 

L is a vector of influence coefficients 
for launch mass perturbations 

A trajectory terminates when a stopping parameter, S, is satisfied. The 
influence of the controls and the launch mass on the stopping parameter is 
represented by the equation 

ois ~ Kldfc 4 r4<S rr,' (2) 
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It is not efficient to interpolate to the correct value of the stopping 
parameter when perturbing the controls. The reason for this is that each 
attempt at interpolation requires the same computation that is made to 
integrate the motion for the entire stage. Therefore, the sensitivity 
matrices will be evaluated for a fixed value of the mass ratio, $ . They 
will then be converted to sensitivities with respect to a fixed stopping 
parameter. The subscripts $ and $ will be used to indicate whether 
the matrix refers to a fixed mass ratio or fixed stopping parameter con- 
dition. 

The following equations relate terminal perturbations at a fixed stopping 
parameter to the perturbations &l a fixed mass ratio. 


* •‘t’l, * r'l f uf (3) 

as| £ * d.s | j, * s'|j, us (h) 


where 


y' 



o' - 

b " JLS 


and df is the change in the nominal mass ratio required to meet the stopping 
condition in the presence of a control perturbation. 

By definition, J,s|g*0 . Solving Eq. 4 for <*■? , one obtains 


Substitute Eq. 5 into Eq. 3 to obtain 





(6) 


Eqs. 1 and 2 are to be evaluated for fixed mass ratios. Using the proper 
subscripts, substitute then into Eq, 6 to obtain 


cLf g - Kg<L& + LgJLm* 


(T) 


where 


«■ Kp — 


Vp' 

fi' 


L'C " i-o “ 


*V' 

S' 


N, 


It is now desired to find the vector Ae which will produce a desired 
change in the presence of a launch mass change A™ . Further- 

more, a weighted sum of the squares of the changes in the control parame- 
ters should be minimized, i.e., one desires to minimize 

i & Z <L e 

while satisfying Eq. 7- Z is a diagonal matrix of weighting functions. 
Use Lagrange multipliers to form the quantity 


V - At, Z Aa It-M* - L r A "-i> 


( 8 ) 
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Write the first order perturbation equation for V 


gv = (2 ItZ -yOKg) g(Afe) (9) 

S V must be zero for arbitrary changes in Afe . The coefficient of S(JLt) 
must therefore be zero. Solving for Jit gives the following result. 

dfc r £ (10) 

Z 

To determine /a , substitute Eq. 10 into Eq. 7* 

f T r OHs " L S <»-*";) (11) 

Substitute Eq. 11 into Eq. 10 to obtain 

ie = ^"K/(K s E"K s T r(if| £ (12) 

In the program, the following notation has been used 


= ^KZ 

K ff <= ^KS 

L S -- IS)\_Z 

L s - &P-S 

tA s ~ 


Nj, = 

M s - ^NS 

Y' * PSIPR 

$' = SPR 

^Xs' = PDS 

K s £"'z/ * A A 

Z*V(K«Z'V)'= ZKA 

(K,z"Z, t )'' - A* 

z"k^(k s z' Kj T ) 1 L s = ZkM. 

Z - '* 2INV 
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K,, Ly, M, , Hj , 4 1 ', and S' are evaluated in the PERTRB sub- 
routine, Kj and Mj, are determined by perturbing the controls, 
one by one, and observing the effect on the terminal constraints. 

L j> and Ny are determined in the same way with the launch 
mass perturbed. V* and S' are evaluated by perturbing the 
mass ratio of the final stage. 

K ^ j> L 5 , and Klj are computed in the MEQ2 subroutine. AA is 
formed and inverted and substituted back into AA. ZKA and ZKAL are 
then formed. 

ZINV is computed at the start of the main program. The actual computa- 
tion of the change in the control vector (Eq. 12) is made in the CtyNC 
subroutine in the form 

DC^N « £KAlV - ZKAL JL^ (13) 

where DC^Kl is the program symbol for de. 


OPTIMIZATION SUCCESS CRITERION 

The optimization success criterion is identical to that used in PRESTO. 

Only the notation has been changed. Thus, Eq. 1 on page 8-l4a in the 
PRESTO report becomes 

= PP&Y Arr) L -4- ^ PP^J dL^V: W 

+ 1 * • 
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An optimization run is successful if the terminal constraints are suf- 
ficiently close to the desired values and if 


( f + <L T(\ * \ "]> (n^r ♦ A fT> r N 

' + T / CURRftMT RUN “ (IlMT RUN 


( 15 ) 


jc 

where dm r = 2 PPAY • dLY: for the fixed launch weight option 

t j'i “1 » 


and 



C IC 

i y pray. a.'v. 

( i'-i 3 


+ PPaY 1 (dWFUEl/ 3.)| (t- P PAY t )"' 


for the fixed final stage option. 


SELECTION OF MASS IMPROVEMENT 


The selection of mass improvement differs from PRESTO in that the improve- 
ment to be asked for is recomputed after every successful run using the 
re lation 


d nn 


WOM 



(l6) 


dm„ 0M is then added to dm^ to obtain d'f, , which is the actual 
change in final mass requested. If a run is unsuccessful, dm MeM is 
cut in half. In the output, DMASD still indicates the first computation 
of made after a successful guidance run. 
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6. EQUATIONS OF MOTION 


AN APPROXIMATE SOLUTION TO THE THREE-DIMENSIONAL EQUATIONS 
OF MOTION OUTSIDE OF THE ATMOSPHERE 


If motion is restricted to the vicinity of a great circle defined by 
the initial velocity vector, it is reasonable to represent the Earth 
as a cylinder. On a cylindrical Earth, motion in the longitudinal and 
lateral directions is uncoupled. Longitudinal motion may be described 
in the usual polar coordinate system vhlle the lateral motion takes 
place under the assumption that the Earth is flat. 

The coordinate system used in PRESTO II is shown on the next page. 

The initial velocity vector lies in the 1 plane. The axis 

is along the local horizontal and the T axis lies along the radius 

y ^ , 

vector, u, v, and v are the components of velocity along the i , i , 

x j 

and I axes, respectively, is the in-plane range angle. A is the 
out-of -plane distance divided by the radius. 

The thrust orientation angles are illustrated on the next page. 9 is 

the angle between the i axis and the component of thrust in the i - i 

X y 

plane. X is the angle between the 7^ axis and the component of 
thrust in the i^ - 7^ plane. ^ is the angle between the thrust vec- 
tor and the local horizontal. It is related to the two control angles 
by the Equation: 
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DIAGRAMS OF COORDINATE SYSTEM 
AND THRUST ORIENTATION ANGLES 
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£ivrv 7J « jtivn* 6 


( 1 ) 


With lift and drag neglected, the equations of motion are: 


a 


X 

m 


MH,7/COn,X Ajrv f 


UV 

r 


( 2 ) 


V = X 

(3) 

w = ~ Yj iLrv X ^ 

M 


(5) 

V 

(6) 

* w 
A= — 

(?) 

m= - i 

(8) 


Sgn 4 is normally +1. It is set to -1 for retroburns because 
0 and X are limited to ±90 degrees. 

It will be convenient to use the mass ratio, instead of time, as the 
independent variable. The mass ratio is defined as: 

J * <9) 

where IPH^ is the mass at the beginning of a stage. The time deriva- 
tive of i is: 


6-3 


dj _ _x_ 

(10) 

With the help of Eq. 10, the equations of motion are converted to: 

(u) 


cL LL 

c C0&X iff t £ m L c 

dJ 

j> + Tr 

JLv 

c airfi 7f f m i c tt* 

cL3 ~ 

jp T r 

olw 

C- C0Tt Aon. X 

dJ ~ 

$ 

<L/S 

c. a 


T r 

JLr 

w:c. 

cLf 

T V 

oLA 

™-c w 

dj> 

' T r 

<L rr\ _ 


df 



U.V 


where c is the exhaust velocity. 


( 12 ) 

( 13 ) 

(14) 

( 15 ) 

(16) 

( 17 ) 


The two thrust orientation angles, & and X > are to vary in accord- 
ance with the "linear -tangent" law; i.e., 


tdM 0 


+ bt 


(18) 
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tam, X - e. + -Pt 


( 19 ) 


a, b, e, and f are constants. 

Using the relation 

w «=• mi - t 


one can write 

ton. e » A + B'i 

U*. X ‘ E Fi 


where A = cu -t- 

T 


E = e + 




B 


be m; 
T 


F 


few; 

T 


With the aid of the following three triangles, the trigonometric 
in the equation of motion cm be written as: 


1 


(Km 7} (XhL y. 


~fl l (A+ B J ) 1 + (E+FJ) 2 -’ 



<KXL7{ Aon.y - 


A +B3 

~fi *(A + B 5)*- + (E+-Fj > ) i ’ 

E + F j 

*j/i + 


( 20 ) 

( 21 ) 

( 22 ) 

terms 

( 23 ) 

( 2*0 

( 25 ) 


6-5 


y i +(e * Fi)*- ' 


The equation for the horizontal component of velocity (Eq. ll) is 
first integrated with the term containing uV neglected. One has, 


" u i 


- - c 



AA 

s yiv^vB j) 2 - 


After integration, one obtains 



where 


Ifl +(a + b^) 1 yTT7 r 7i r ’ + t 4- a(a+b$,) 4-e(e+ Fl t ) 

i 1 + (K l &)* +(& + F f-jj 1 + ^ + E> -*- i + /\(A*B)* E(E + F) 


For 


/» 


A -A 


one has 



{ 


T r 



i) 




Y* 

i 


(i tJL 

IJ 


oLu. 

u 



J.1 J.f 


<Ll 


The double integral is evaluated with the help of the relation 


i 1 1 


( 27 ) 


Using Eq. 27, one obtains: 


^ w ic 

"A = TV~ 


+ *i(VO + ^=T Y 


( 28 ) 
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where 


Y - U 


ab»ef 


If B and F are both zero, the last term in Eq. 28 is: 
c(?,-l)^f 


In order to integrate the equation for vertical velocity, an integrable 
expression for u5 must be found. The approach used here is to split 
the logarithmic term for a into two parts and then to expand the 
second term in a series in 3 . One can write 


(29) 


a * u- 1- 


/immF 



fi+kfaeS)4,Efct-F3) 

^i+A^MS)* E(E + F) +lft+ /V* + E*’ i '|l4(A*.B) 1 -4 (E ♦F) 4 -' 


Note that if the slopes of the control angles are equal to zero, the 
second In term is equal to zero. Furthermore, when B and F are 
not zero, numerical checks for reasonable control angles indicate that 
the term in the round brackets does not move very far from unity. Write 
the term in the round brackets as N/D, where N and D are numerator 
and denominator, respectively. Then: 
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assuming N/D Is close to unity. Using the approximation 
fL + <L' ~ i + —■ 

one obtains for Z the equation: 

Z « (1 +K !♦ + (AB+ E F)(f - 1)* F l ) (i - lj ( 30 ) 


U, can then he written as 


„ _ „ * S ^T lI I 

1 fi TW7V? J/ *' 


X z 
i * d 


(3l) 


After collecting terms, one has 


UL 


= A t + t\ z JLl + 


(32) 


where 


A, = 


__ ^ ^s) 

^4 


A_ = 


A„ = 


^4 


A,. = 


<^A 


and 


a, = 


<L 

| i /a^e* 
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°2 


= (i yrrs r ri T )(A&^E.F) 

Vi+ a l -+ e 1, ' (b^+f*-) 

3 = a 


a^ = 1 4-A(A + Bj •»■ E^E + Fj +^l4 A^+ Yl+-(A V + B) 4 '<-(E + 



1b now readily integrable. 


Eq. 12 becomes: 


cl V 

df 


c(A^bS) 

■5]/T7f7&Jf + (E>F3) tl 


4- 


T 


+ B t + b z 3 +- B 3 'S i + b + , S 8 + e> s i 




(33) 


where 

B x -KA‘ 

B 3 = K(A 5 t ^A 1 A 4 ) 


= Z K A 3 A 4 
= K A/ 

Bg =r ZkA i A i 


k 


m;c 

Tr 


B^ = K Aa, 

B g = ZKA*A 3 

b 9 = zka x a 4 


Integration of Eq. 33 yields: 

e E3> - aa/tv £ m: c 


/ \ 



where 

* e «[r (J*% ~i) 4 "£] + -•V) 4 - t 


Using Eq. 2J once again, the equation for f* is integrated to give: 



If B and F are zero, the term containing Y becomes: 
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In a similar manner, the equations for v and A are integrated to 
give: 


- w t 


c E ^ 

7 L *• A * + E*-' 


cF ^ Y 
7 B + F T 


(36) 


, A m;c r, / e B (BE - ^ 

A - Ai = W L j + + p *-) TJi T (37) 

~~ (j 1 + (A+ B$l Y + (E + F ^ ) t ~- 'yi + CA + B^E + F) 1 -)* w £ (.V 


+ A 
&* 

If B and P are aero, the term containing Y becomes: 


Y 1 -v A r + E 1 "' 


It will be recalled that the term containing uv in the differential 
equation for horizontal velocity was neglected. The influence of this 
term is accounted for in an approximate manner by adding to ^ the 
average value of — uv multiplied by the interval of inte- 

gration. This additional term is: 


ThiB approximate solution to the equations of motion is repeated for 
each powered stage with the final conditions at the end of the first 
stage becoming the initial conditions for the second stage. If a 
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coast precedes a powered stage, the axis of the cylinder Is rotated 
so that there Is no lateral notion at the beginning of the powered 
state; i.e., A- t and w; are set to zero. 

The value of r that is used in the solutions is obtained by averag> 
ing the initial and final radius for a stage on the last iteration, 
g. is computed using this average value of r . 
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CONVERSION FROM CYLINDRICAL TO GEOGRAPHIC COORDINATES 


At the end of each powered stage, it is necessary to determine the posi- 
tion and velocity of the vehicle with respect to the spherical Earth. 
This computation is done in the OONVRT subroutine. 

At the beginning of the first control segment and after each coast we 
have available the inertial azimuth, , the latitude, \ , and the 
longitude, 'T' . At the end of each powered stage, we know u , v , 
w, fi, and A. It is necessary to compute inertial velocity, 

V T , inertial flight path angle, T r , inertial azimuth, latitude 
and longitude. 

Except for the symbols mentioned above, the terminology used in this 
subroutine is not related to the rest of the program. 

Inertial velocity is easily found from the Equation: 

V r = * v 2 +• w 2 ’ (l) 


Inertial flight path angle is found from: 


( 2 ) 
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Spherical trigonometry is required to determine the remaining variables. 
Consider the following triangle: 


North Pols 



Poimt 




The computation is broken into two parts. We first determine the latitude 
and longitude at point A which is on the great circle defined by 
the initial velocity vector. 7 takes on one of two values depending on 
at the initial point: 

If % + , 7 - 4 *! 

if > tt , 7 = £ir - ‘Vi 

S -E-X 


Using spherical trigonometry, one obtains: 

C&O, CK - Cx>u/ £ CasC,' ft <*1™ £ 7f 


(i) 
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&vrv oc 

A^n, 


Acre ST 


Compute Aun. oc 


from (l) as: 


AUsn- oc - fi - (UKls^ct ’ 


Use Eqe. (1) and (3) to compute: 


ot 



If tauYi oL > 0 , 

If XOSTL CC <C 0 f 
Knowing cC , solve for 
Aim- St = — 


0<L | 

occ TT 

Z 

AOrc $ T 

aUi Acm^/S 
din oc 


Then, Cod/ ST = ~l 1 - 


£Tis always assumed to be less than 


The latitude at point A, X A , is: 


X 


ir 


K ~ T 


- °C 


and the longitude is: 


( 7 ) 




+ if ^ 'Tr 
- if Vj7ir 


Knowing ST" , one can solve for c using 

Co& C - — Ccrti/ ST Co&Tf 4- xiwt/ ST ^n. ^ CCrrV S 
= y j - Co& Z C ' 

and 


( 8 ) 



c « 

, -i / a^rt-c \ 

(9) 

If 

e ^ <3 

7T 


If 

i«n.c ^ O 

tr 


Then 

tr 

ii 

=1 

i 

o 

(10 


The inertial azimuth at the end of a stage also depends upon the out- 
of -plane velocity, w . Let 

d = (ll) 
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The inertial azimuth is then found from! 


+<1 1$ ^ ^TT 

- 2/TT - + cl I-f Y- ^Tr 

*A 


( 12 ) 


The second part of the calculation is similar to the first except that 
TV replaces . 



Define an angle m as follows: 

m « ^ ^ Afn d if ^ -tr 

m = ^'Tr - ‘Vj + Ayrv<i' if ‘7'jr 

If m >■ ^TT, subtract 2 "TT from it . 
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Define 7? 


as follows: 


Tj = m 


if 


7 

5 

and 

CtfU/ oC 
&un. ct 


- z'Tr-m 


IT \ 
T " a A 


if tn 7ir 


C6M 5 «a,A 4- xiArv 5 dinft/v C^O/7^ 
jjl - CM z oc' 


ai = 


„l( AttVcL 
^ a/yL \ CM/ oL 


) 


(13) 


Use the quadrant check that follows Eq. 4: 
Then 

. „ UX/rv 7? dim, A. 

^-f r = — ^ — 


CM/ ST - ^ i ~m*i z ST" ' 

and 

c _ . -x / xtAn- ST - ^ 

° ' = ' z *~' £-r ; 


The latitude at the end of the stage is: 




and the longitude is: 

« Ti - OJAt 1 ST 


+ if m ^ IT 


- if m ^ tt 

where At is the time from the initial point. 


(14) 


(15) 


(16) 
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CONVERSION OF PRESTO II CONTROL ANGLES TO PRESTO CONTROL ANGLES 


At the end of the PRESTO II iteration, we have two control angles, 0,, 
and Xg> defined with respect to the PRESTO II coordinate system. It is 
necessary to convert these angles to the PRESTO control angles, T| and X, 
for use in generating the initial integrated trajectory. 

The two coordinate systems are related as shown in the following diagram. 
The subscript 2 refers to the PRESTO II system, y is the flight path angle. 



6i|r is the angle between the great circle path defined by the initial velocity 
vector in the PRESTO II system and the horizontal component of velocity in 
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4 


the PRESTO system. The lateral motion in the PRESTO II system, described 

by the angle A Is assumed to be small. This leads to the approximation 

that the i and i axes are parallel when y 1 b zero, 
z y 2 

The orientation of the thrust vector in the two coordinate systems 
is shown below. 


i 

z 

4 


i 


y 



PRESTO 




PRESTO II 


The thrust vector in the PRESTO II system is resolved into components 
to give 


T = T i T +i T + i T 

“ “ y 2 y 2 z 2 ^2 


i 





L 

*2 

Xg 

T 

x 2 

= 

COS 


cos 

*2 

T 

y 2 

ss 

sin 




T 

z 2 

ss 

cos 


sin 

Xg 

tan 


= tan 

0g cos 


(i) 


The components of thrust in the PRESTO II system, resolved along the 
PRESTO coordinates, is given by 
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T = T | i x |_ _T X ^ sin 6 i/ + T^ cos 6+ J + cos cos y + sin y 


+ T sin 6* cos y + i -T cos 6* sin v + T cos y 
z 2 'J z l *2 ?2 


- T z sin 6t|f sin yj j- 


The components of thrust in the PRESTO system, resolved along the PRESTO 
coordinates, is given by 


T « T j^i x cos T] sin x + i cos T| cos x + i 2 sin T)J 


The coefficients of the unit vectors in Eqs. 2 and 3 must be equal. One 
obtains 


sin T) = -T cos sin y + T cos v - T sin sin y (4) 

*2 y2 z 2 


cos = / 1 - sin T) 


• 1 r sin T) 
V cos T) 


sin x = 


cos T) L Xg 


j -T sin 6* + T cos 
L Xg z 2 


cos X = y"l - sin 5 X 


■ 1 f' sin x 
V cos x 


These calculations are made in the PCAL subroutine. 


Hie angle 6t|r is computed as 


« t - t, 
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where i is the local azimuth and t R is the instantaneous azimuth along the 
great circle defined by the initial inertial velocity vector. is computed 
after every step of integration using the relation 



Vi * *» (DT> 


where 


r R 


Vj cos Vj sin sin X 
r cos X 


( 11 ) 


NAS A- Langley , 1967 CR-686 
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conducted so as to contribute ... to the expansion of human knowl- 
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shall provide for the widest practicable and appropriate dissemination 
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